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Abstract. We consider random non-normal matrices constructed by remov¬ 
ing one row and column from samples from Dyson’s circular ensembles or 
samples from the classical compact groups. We develop sparse matrix models 
whose spectral measures match these ensembles. This allows us to compute the 
joint law of the eigenvalues, which have a natural interpretation as resonances 
for open quantum systems or as electrostatic charges located in a dielectric 
medium. 

Our methods allow us to consider all values of /3 > 0, not merely /3 = 1, 2,4. 


1. Introduction 

The main objects of investigation of this paper are Dyson’s circular ensembles of 
unitary random matrices, as well as the orthogonal and compact symplectic groups 
equipped with the Haar measure. To give the flavor of the results, let us restrict 
our attention in this introduction to Dyson’s circular ensembles only. 

In 1962 Dyson introduced three ensembles of unitary random matrices {CUE, 
COE, and CSE) to model complex physical systems (e.g., evolution operators of 
closed quantum systems) corresponding to various physical symmetry classes. We 
begin with the definitions. Note that the CSE ensemble is best explained through 
the use of (real and complex) quaternions. We review this material and introduce 
the relevant notations in Appendix [A| 

Definition 1.1. Dyson’s three circular ensembles are: 

(CUE) The circular unitary ensemble CUE{n) is the set {group) of all nxn unitary 
matrices U(n) endowed with the Haar measure {which comes from the group 
structure o/U(n)). 

(COE) The circular orthogonal ensemble COE{n) is the set of all nxn symmetric 
unitary matrices with the measure induced by the Haar measure on U(n) 
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via the mapping 

U(n) ^ COE{n) 

U hA U'^U. 

(CSE) The circular symplectic ensemble CSE[n) is the set of all n x n complex 
quaternionic matrices that are both unitary and self-dual. The measure on 
this set is taken to be the one induced by the Haar measure on l[J(2n) via 
the mapping 

U(2n) ^ CSE{n) 

U ^ <L-\U^U). 

As argued by Dyson, COE is relevant in most practical circumstances, with the 
two exceptions of systems without time-reversal invariance (when CUE should be 
used) and odd-spin systems with time-reversal but without rotational symmetry 
(when CSE should be used). 

Dyson computed that the eigenvalues of these ensembles are jointly distributed 
on 9D" := {e*® : 0 < 0 < 27r}" proportionally to 

oc d0i...d0„, 

j<k 

where /3 = 1,2,4 for COE{n), CUE(n), CSE{n), respectively (for CSE{n) each 
of the n eigenvalues is of multiplicity 2). 

For each of the listed ensembles we consider their truncations: choose at ran¬ 
dom a matrix from the ensemble and delete any row and the corresponding column 
(quarternionic row/column in the CSE case). These truncations appear naturally 
as the evolution operators of open quantum systems, and the eigenvalues of these 
truncations can be thought of as the scattering resonances of a random quantum 
system attached to a “lead”. This and other physical applications of truncations 
have recently gained a significant popularity in the mathematical and physical lit¬ 
erature. We refer the interested reader to excellent review papers na [ 111112 ] and 
references therein. 

The unitary truncations have connections in combinatorics (see m) and to the 
theory random Gaussian analytic functions (see HIIIIIHI)- 

To be able to systematically analyze the truncations, it is natural to start with 
finding the joint eigenvalue distribution. For the truncations of CUE{n 1) it has 
been computed in the physics literature by Zyczkowski-Sommers [35] (see also [291 
HHETj). The result is that the eigenvalues are distributed in D” = {z S C : |z| < 
1}" according to 

n 

j,k=l j<k 

where d^ is the two-dimensional Lebesgue measure on D. In this paper we solve 
this problem for COE and CSE. 

Theorem 1.2. The eigenvalues zi, ..., of COE{n-\-l), CUE{n-\-\), CSE{n-\-\) 
ensembles with one row and column removed are distributed in D” according to 

n 

n - ZjZk)^~^ Y[\zk - Zj\^ d^zi.. .d^z, 

j,k=l j<k 


( 1 . 1 ) 
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(a) (b) (c) 

Figure 1. A realization of random eigenvalues for truncations of: 
(a) U(n+1) with n = 301; (b) 0(n+l) with n = 301; (c) USp{n+l) 
with n = 151. 


with /3 = 1, 2,4 respectively. 


Remarks. 1. For CSE each of the eigenvalues is of multiplicity 2. 

2. In fact, we show that for any 0 < /3 < oo ([nj is the eigenvalue distribution of 
the circular /3-ensemble of Killip-Nenciu (see Remark 1 after Proposition 4.2) 
with the first row and column removed. 

3. /3 = 2 is quite special here as the process becomes determinantal [55] . 


One of the ingredients of the proof is to reduce these truncations to a special 
sparse matrix form (namely, CMV form, see Dehnition 2.3) that depends on only n 


independent coefficients with explicit distributions. This is done in Proposition |5.2| 
and could be of interest on its own for simulation purposes. 

Another way to look at this is that we show that the eigenvalues of the trun¬ 
cations coincide with the zeros of orthogonal polynomials on the unit circle with 
random recursion coefficients that are independently distributed according to the 
explicit distributions (5.1). 


Similarly, we develop matrix models for truncations of the orthogonal group 
(and, more generally, of the real orthogonal /3-ensemble) and of the compact sym- 
plectic group. We establish connection to zeros of orthogonal polynomials and 
compute the distribution of the (independent) recurrence coefficients. Finally, we 
compute the distribution of eigenvalues of these ensembles (equivalently, distribu¬ 
tion of zeros of these random orthogonal polynomials), with the exception of the 
eigenvalue distribution of the truncated compact symplectic group, which remains 
an open problem as of now. For a reader who is wondering why the compact 
symplectic group seems special among all the other classical ensembles, we remark 
that deleting one quaternionic row and column corresponds to the deletion of two 
complex rows and columns. Truncation of more than 1 row and column can also 
in principle be attacked with our methods, which leads to a CMV model with 
matrix-valued Verblunsky coefficients. However, the Jacobian computation is not 
technically manageable to us at this point. We leave this as an interesting open 
problem. 

The joint eigenvalue distribution for truncated 0(2fc) was found earlier via ad 
hoc methods by Khoruzhenko-Sommers-Zyczkowski in [20] . Since the completion 






4 


TRUNCATIONS OF UNITARY ENSEMBLES 


of this work, Forrester [9] reported a formula for the joint eigenvalue density of 
truncated US'p(n) ensemble. 

Our methods can be applied to random ensembles with the so-called non-ideal 
coupling (Fyodorov-Sommers [TS], Fyodorov [T^], and Fyodorov-Khoruzhenko m)- 
We produce /3-matrix models and compute the joint eigenvalue density. 

After establishing the joint distribution, it is natural to ask what can be said 
about the correlation functions and their large n behavior. In our upcoming pa¬ 
per m we investigate the microscopic limit of the 1-point correlation function of 
the circular /3-ensemble (0 < /3 < oo). In fact, we show that the same limit is 
obtained for a wide class of subunitary random matrices. 

The organization of the paper is as follows. In Section we introduce the rest 
of our unitary random matrix ensembles (the compact groups) and define CMV 
matrices and the CMV-fication algorithm. In the rather lengthy Section we 
compute the random spectral measure of each of the ensembles we study. Their 
eigenvalue distributions are well-known of course, but we also need distributions of 
the eigenweights of the spectral measures. In Sectionj^we present the matrix models 
that are obtained by CMV-fying each of the (non-truncated) unitary ensembles. 
This is originally the idea of Killip-Nenciu [52], and we extend this to all the other 
unitary ensembles. 

In Section [5] we derive CMV matrix models for each of the truncated ensembles. 
In Section [^ we compute the joint eigenvalue distribution of the truncations. In 
Section we discuss matrix models and compute the eigenvalue distribution for the 
non-ideally coupled ensembles. 

In Section]^ we show that the eigenvalues of the truncated Dyson’s (and, more 
generally, circular /3-) ensembles admit a log-gas interpretation as charged particles 
located in a dielectric medium. 

In Section]^ we introduce a symmetric variant of the CMV model, which allows 
us to find a canonical system of representatives of each conjugacy class that actually 
lie in the sets COE and CSE; the traditional CMV matrix representation has no 
such symmetry properties. Finally, in Appendices we collect the information about 
quaternions, basics of the theory of orthogonal polynomials on the unit circle, and 
Jacobian determinants. 


2. Preliminaries 

2.1. Compact groups and Haar measure. As discussed in the Introduction, 
apart from Dyson’s circular ensembles defined in Definition o we also study the 
classical compact groups. For the notation related to quaternions and quaternionic 
matrices, see Appendix [A] 

Definition 2.1. The unitary, orthogonal, special orthogonal, and compact symplec- 
tic groups are defined, respectively, by 

U(n) = {Q G : Q^Q = QQ^ = /„} , 

0(n) = {Q G : Q^Q = QQ^ = /„} , 

SO(n) = {Q G : Q^Q = QQ^ = /„,detQ = l} , 

USp{n) = {Q G : €{Q)^Q:{Q) = <t{Q)€{Q)^ = hu} ■ 

Here is the nxn identity matrix; 1 and ^ stands for the Hermitian conjugation 
and transposition, respectively. 
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By the general theory, each of these groups possesses the Haar measure, which 
is the unique (normalized by 1) measure that is invariant under the group action. 
We will use the same notation l[J(n),0(n),SO(n),US'p(n) for the group itself and 
for the random ensemble. 

Next we describe the joint law of the columns of matrices from these groups. 
This is well known, as is the method of proof we employ. We give details only in 
the US'p case where the argument contains an additional subtlety. 


Proposition 2.2. (a) The first column vi of a matrix from 0(n) is distributed 
uniformly on the sphere {a; € M" : ||a;|| = 1}. For 2 < k < n, the k-th column 
Vk is distributed uniformly on the subset of this unit sphere that is orthogonal to 

V \, . . . , V}^ — \ . 

(b) The first column vi of a matrix from U(n) is distributed uniformly on the 
sphere {z G C” : ||z|| = 1}. For 2 < k < n, the k-th column Vk is distributed 
uniformly on the subset of this unit sphere that is orthogonal to vi,... ,Vk-i- 

(c) Let Q G USp(n). The first column vi of (L{Q) is distributed uniformly on the 

sphere {z € : || 2 || = 1}. The second column is Zvi, where Z is (A.3). For 

2 < k < n, the {2k — \)-th column V 2 k-i is distributed uniformly on the subset of 
this unit sphere that is orthogonal to vi,..., V 2 k- 2 , and V 2 k = Zv 2 k-i- 

Proof. As noted above, we give details only for the case of US'p(n). 

Form V with columns ui,..., V 2 n described by the procedure in the Proposition, 
and let us show that V is indeed £(Q) for a Q S l]Sp{n). By the construction, 
ui,... ,V 2 n form an orthonormal basis of C^" since v _L Zv holds for any vector 
V. Thus V is unitary. Moreover, V consists of blocks each of which has the 
real quaternionic form, see (A.2). Therefore C~^{V) belongs to the group (not 


ensemble yet) USp{n). Let us show that for any non-random W chosen from the 
group £(U5'p(n)), the distribution of V and WV are identical. This will prove that 
€~^{V) was generated according to the Haar measure by the uniqueness of the 
measure invariant with respect to the group action on l]Sp{n). 

Denote the columns of WV by ui,... ,U 2 n- First note that ui = Wvi is a 
vector uniformly distributed on the unit sphere {z G : ||z|| = 1} since Vi is. 


The second column U 2 = Wv 2 = WZvi = ZWvi = Zui (we used (A.4)), which 


agrees with our construction. Now inductively, consider U2k-ii which is orthogonal 
to ui,..., U 2 k- 2 - Consider any 2n x 2n unitary U that maps spanjui,..., U 2 k- 2 } 
to itself. Then W~^UW maps spanjui,... ,V 2 k- 2 } to itself. By the construction 
this means that V 2 k-i and W~^UWv 2 k-i are distributed identically. Thus U 2 k-i = 
Wv 2 k-i and Uu 2 k-i = UWv 2 k-i are distributed identically. Since U was arbitrary, 
we conclude that U 2 fc-i was distributed uniformly on the subspace of C^" orthogonal 
to Ui, . . .,U2k-2- 

Finally, U 2 k = Zu 2 k-i as before, and this completes the induction. □ 


2.2. CMV matrices. Let us now give the definition of CMV matrices. They 
are named after Cantero-Moral-Velasquez [1] (also discovered earlier by Bunse- 
Gerstner-Elsner M)- For a comprehensive discussion, see [5^ . 


Definition 2.3. Given coefficients oq) • ■ • ,0(n-i with Ofe £ D for 0 < k < n — 2 
and Un-i G D, let pk = -\/l — and define 


— 


ttfc Pk 
Pk —Oik 


( 2 . 1 ) 
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for 0 < k < n — 1, while S_i = [1] and E!„_i = From these, form the n x n 

block-diagonal matrices 

/: = diag(So,S 2 ,S 4 ,...), = diag (S_i, Si, S 3 ,...) • ( 2 . 2 ) 

The matrix 


• • ■ : 1 ) ■— 

is called the CMV matrix corresponding to ao,... 


Remark. Note that C(ao,..., Oin-i) is unitary if and only if a„_i € 9D. If an-i € 
D, then it can be shown that C(ao,... , 0 ^- 1 ) is a subunitary operator (C^C < 1) 
with defect operator of order 1 (see [33] as well as my Usually the name “CMV 
matrix” is reserved for the cases when a„_i S (9D, and if an-i G D then these 
matrices are referred to as the “cut-off CMV matrices”. In this paper we will use 
the term “CMV matrix” for all a„_i G D. 


Let us now define what we will refer to here as the CMV-fication algorithm. 
Suppose we are given an n x n unitary matrix U (with complex entries). If the 
vector Cl := [1,0,0,..., 0]^ is cyclic for U then (see [31IH|32|) applying the Gram- 
Schmidt orthonormalization procedure to ei, Uei, U~^ei, C/^ei, 17“^ei,... produces 
a basis in which U has the CMV form C with a„_i G 9D. 

Similar construction works if U is instead a unitary operator on with 

the cyclic vector ei. Then one ends up with an infinite CMV matrix C(ao, oi,...) 
with aj G D for all j > 0. This matrix is defined as CM, where C and M are 


semi-infinite block-diagonal matrices ( 2 . 21 , where are as in ( 2 . 1 ) for all fc > 1 . 

Recall that by the spectral theorem, if U is a unitary operator on or 

on ^^(Z+) with a cyclic vector ei, then its spectral measure is defined to be the 
probability measure on 9D uniquely determined by 


(ei,U’ 


^ei) = / 
^0 


Am9 


dp.{0), 


(2.3) 


where {v,u) = vA. In fact, we have an isometry i of Hilbert spaces L^{ix) and 
£^(C") or £^(Z_|_) defined by 

i hA U”‘ei, mGZ. (2.4) 


Now note that during the Gram-Schmidt procedure above, ei was not changed, 
and therefore our CMV-fication algorithm can be described as 

u = WCW^ 


for some unitary W with lUei = W^ei = ei. 


(2.5) 

This means that not only the 


eigenvalues of U and C coincide, but also that the spectral measures (2.3) of U and 
C with respect to ei are identical. In fact, two unitary CMV matrices with identical 
spectral measures necessarily coincide, see |23l Prop 3.3]. 

Thus from the physical point of view, one may want to think that given a physical 
system modeled by a unitary matrix U with a chosen vector Ci (viewed as the initial 
state of the system) there is a unique choice of basis in which the matrix attains 
the “minimal complexity” CMV form. 

Note however that in the case of systems with time-reversal invariance {COE and 
CSE) the associated CMV form will not possess the same symmetry properties. In 
Section we introduce a variant of the CMV form that does possess the requisite 
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symmetry. Moreover, the matrix used to conjugate to this new canonical form 
belongs to the othogonal/symplectic group. 

From the mathematical point of view, one of the advantages of reducing a uni¬ 
tary matrix or operator U to the CMV form, apart from the latter being sparse and 
having the same the spectral measure, is that there is a direct connection to the 
theory of orthogonal polynomials on the unit circle. Indeed, characteristic polyno¬ 
mials of C{ao, ..., Q!fc_i), k = 1,2,..., coincide with the orthogonal polynomials 
associated to the spectral measure fi of U with respect to Ci (see Appendix [B|) : 

det{zlk - C{ao, ■■■, Ok-i)) = ^k{z). (2.6) 


In particular, this means that these polynomials satisfy Szegd’s recurrence (B.l). 
We stress that the recurrence coefficients in (B.l I are precisely the coefficients from 
the CMV matrix. We will refer to them as the Verblunsky coefficients. 

We collect all the needed facts from the theory of orthogonal polynomials in 
Appendix]^ For more details the reader is referred to the monographs HE]- 


2.3. Block CMV matrices and matrix-valued spectral measures. For the 

quaternionic ensembles CSE{n) and l]Sp{n) we will have to deal with the 2x2 
matrix-valued spectral measures and with block CMV forms which we define next. 

Recall that p. is an / x / matrix-valued probability measure on 9D if it is a 
countably additive mapping from Borel subsets of 9D to positive semi-definite I x I 
matrices, with the normalization condition /r(9D) = The I x I spectral measure 
of an In x In unitary matrix U with respect to vectors ei,e2,...,e; is the I x I 
matrix-valued measure dp with entries dpjk (1 ^ j,k < 1) determined by 

(e„C/”^efc)=/ E^^dp,k{0). (2.7) 

Jo 


Block CMV matrices can be defined as before, by replacing the scalar Verblunsky 
coefficients aj G V) hy I x I matrices satisfying a^a < Ii and pj = {Ii — 

Under the natural cyclicity condition, any unitary nl x nl matrix U is unitarily 
equivalent to a block CMV matrix W^UW = C(ao, ■ ■ ■, On-i) with the last co¬ 
efficient satisfying al^_ian-i = Iixi- Moreover, the I x I matrix-valued spectral 
measure (2.7) with respect to ei,e 2 ,...,e/ is preserved (which is to say, W fixes 

ei,...ei). 

Now, given a matrix Q from CSE{n) or l[J5'p(n), by its spectral measure we will 
call the 2 x 2 matrix-valued spectral measure of the 2 n x 2n complex matrix €.{Q) 
with respect to ei,e 2 . The matrix £(Q) can be reduced (with probability 1) to 
C(q;o, • ■ •, Q!„_i) with 2x2 Verblunsky coefficients aj. The 2x2 spectral measures 
of £(Q) and C with respect to ei, 62 coincide. 


3. Spectral Measures of Unitary Ensembles 

In this section we will derive the distribution of the random spectral measures 
of each of the ensembles defined above. 

Let us denote the Vandermonde determinant by 

A{zi,...,Zn) = {Zk-Zj). 

l<j<k<n 
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3.1. Dyson’s circular ensembles. 


Proposition 3.1. The spectral measure of COE{n) and of CUE{n) is 


n 



E/ ’ 

7 = 1 

(3.1) 

where 



0 < < 271 

(1 < j < n), 

(3.2) 

n 

E/ Tj — 

i=i 

Tj >0 (1 < j < n), 

(3.3) 

and the joint distribution of e*®’ ,.. 

. ,e*®",/xi, .. .,pn-i is 


z„,f3 1 ^ 271 

■^''^11^7 d/3,...dpr.-l 

(3.4) 


i=i 


with /3 = 1 for COE and (3 = 2 for CUE. Here Zn,/ 3 , Z'^ ^ are the normalization 
constants expressed in terms of the gamma function as 

^ _ r(i/3n + i) _[rmr 

[r(i/3 + i)]”’ n\M ■ 

The 2x2 spectral measure of CSE{n) is 


E 

i=i 


fij 0 
0 


(3.5) 


where p,j and 9j follow the law (3.4) with 13 = A on the domain (3.2)-(3.3). 


Remark. In fact, the probability distribution of the eigenvalues in (3.4) with any 


0 < /3 < oo has a natural interpretation as the Gibbs distribution for the classical 
Coulomb gas on the circle at the inverse temperature (3. This distribution for 
any 0 < ,5 < oo was realized as the eigenvalue distribution of the family of random 
matrix models called circular /3-ensembles in Killip-Nenciu [52]. See Proposition |4.2| 
below. 


Proof. The distribution of the eigenvalues is well-known (see Dyson’s original pa¬ 
per (Sj). 

Distribution of the eigenweights for /3 = 2 was computed in |22| Proposition 3.1]. 

Let us compute the distributions of the eigenweights for /3 = 1. Suppose U G 
COE{n), and let its eigenvalues be given. It can diagonalized by 

a real orthogonal matrix U = ODO^ (see (5] Thm 4]; or, more generally, [T71 
Thm 4.4.7]). The matrix O of eigenvectors is unique up to a multiplication of it 
by a diagonal matrix S with ±1 on its diagonal. If we choose each of these signs 
by random independently with equal probability, then the map U i— OS becomes 
well defined, and OS becomes a Haar distributed orthogonal matrix. Indeed, the 
invariance property of COE{n) ([221 Thm 10.1.1]) implies that the distribution 
of OS and WOS are identical for any orthogonal IT, and now we can apply the 
uniqueness of the Haar measure on 0(n). Now Proposition [2.2[ a) says that the 
first row X = {xi,... ,Xn) of OS is uniformly distributed on the unit sphere S^~^ — 
{a; G M” : I 
the proof. 


= 1}. But by (2.3) and (3.1), g,j = Now Lemma C.l a) finishes 













TRUNCATIONS OF UNITARY ENSEMBLES 


9 


A similar approach works for /3 = 4. Let U G CSE{n). Suppose the 2n eigen¬ 
values of €{U) are i e*®", e*®" (each repeated twice). By [SJ Thm 3] 

U can be reduced to the diagonal form U = BDB^ by a symplectic matrix 
B G USp(n), where €(-D) is the 2n x 2n complex diagonal matrix with the di¬ 
agonal ,..., e*®", e*®"). B is defined uniquely up to multiplication by a 

diagonal matrix S with unimodular real quaternions on its diagonal. If we choose 
those by random independently and uniformly on l[JS'p(2), then by the invariance 
property of CSE{n) (see [5SJ Thm 10.2.1]), BS becomes Haar distributed random 
matrix from USp{n). 

Denote the first quaternionic column of {BS)^ by [qi,..., with 




aj Pj 

Pi ^3 


with aj,/3j G C, 1 < j < n. Now C/"* = {BS)D"^{BS)^ and a simple calculation 
implies 


(ei, e:([/)-ei) = ( 62 , e{ure 2 } = + PjlPe 

i=i 

(ei,e:([/)'"e 2 ) = (e 2 ,e:([/)'"ei)= 0 . 


imOi 


\^J \ ^ Ih'j 

is [ai,/3i,. 


Therefore the 2x2 spectral m easu re \2.1\ of U is indeed of the form (3.5) with pj = 
By Proposition 2.2 ]c), the first complex column of €{{BS)^ 
is uniformly distributed on the unit sphere of C^”. 


that 

Then 


by Lemma C.l ]e), pj’s are distributed according to (2n — 1)! n?=i • ■ • dpn-i 


on (3.3) as required. 


□ 


3.2. Orthogonal group. 

Proposition 3.2. (a) The spectral measure o/SO(2n) is 

n 

(3-6) 

where 


0 < 9j < TT 

P<j< 

n). 

(3.7) 

n 

Tj— 
i=i 

Pj > 0 

(1 < j < n), 

(3.8) 

and the joint distribution of 

iOi 

c 5 • ■ • ? 

... ,/x„_i is 



|A(2cos6»i,... ,2cos6'„)|^ d9i .. .d 6 »„ x (n - l)!d^i. ..dpn-i- 
(b) The spectral measure of 0{2n) \ SO(2n) is 

n—1 

h'TiSl + Pn+lS-l, 

i=i 


(3.9) 
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where 


Q < 9j < TT (1 < j < n — 1), 

n+1 

= > 0 (l<j<n + l), 

and the joint distribution of , • ■ ■, , /ii,..., /i„ zs 


(3.10) 

(3.11) 


|A(2cos6'i, ... ,2cos6»„_i)|^ sin^ 6j d0i... d0n-i 


1=1 


(«-!)! 1 
7r 


d/ll ■ • • dun- 


(c) The spectral measure o/SO(2n+ 1) is 

n 

1 = 1 


with {3.7), (3.11); and the joint distribution of , fii,... 


2^ 

nlTT^' 


n 2 

|A(2cos0i,..., 2 cos0„)|^ sin^(^) d0i ... d0„ x -^=L=d/ii. 

,=i r(2) M.+1 

(d) T/ie spectral measure ofO{2n + 1) \ SO(2n + 1) is 


(3.12) 
; i^ 

■dpLn. 


sMl ) + l^n+1^-1) 

1=1 


with (3.7), (3.11), and the joint distribution of e''^^,... , pLi,... 


2” 

nlTT’' 


|A(2cos6»i,...,2cos6»„)|^]Jcos^(^)d6»i...d6»„ x ^ 7 p=^Mi ■ 

j=i ^^ 2 ) 


(3.13) 

, /lyi is 
.. dfiji- 


Proof. The distributions of the eigenvalues is well-known and can be found in, 
e.g., [33] and [55], so we just need to justify the distribution of the eigenweights. 

(a) is proven in [331 Prop 3.4]. 

(b) , (c), and (d) can be proved in the similar fashion. Indeed, one can reduce a 
random matrix to the block diagonal 


diag ^ 

cos 01 
— sin 01 

sin 01 
cos 01 

7 • • ■ 5 

diag ^ 

cos 01 

sin 01 

■ 

— sin 01 

cos 01 

, . . . , 

diag ^ 

cos 01 

sin 01 


— sin 01 

cos 01 

, . . . , 


cos 0n-l 

— sin 0n- 

COS 0n 

— sin 0n 

cos 0n 

— sin 0n 



respectively. 

This shows that the vector of eigenweights (/ii,... ,/in+i) for the cases (b), (c), 
and (d) is distributed as in Lemma C.l (c), (d), (d), respectively, which finishes the 
proof. □ 
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3.3. Compact symplectic unitary group. 

Proposition 3.3. The 2x2 spectral measure ofT}Sp{n) is 

n n 

i=i i=i 


(3.14) 


where 


W, = 


Tj VJW e 


ii/jj 






, wf = 




-y/jrjv]e 

Tj 


itpj 


Here 


0 < < TT (1 < j < n), 

0 < -ipj < 271 (1 < j < n), 

n n 

E + E ^ (1 < j < 

i=i i=i 

and the joint distribution of ..., e*®", , • ■ •, , /ii,..., /i„, ui,..., u„_i zs 


|A(2cos0i,..., 2 cos0„)|^ sin^ 9j d9i... d9n 


i=i 

X ^ X (f^n-^V-dfii. ..dpLndvi. ..dv^-i- 

Proof. The distribution of the eigenvalues is a classical result. As for the matrix 
eigenweights, first write a matrix U G USp(n) as U = BDB^, where £(4?) = 
diag(e®®®, ,..., e*®", e“®®"). Repeating the arguments from the proof of Propo¬ 

sition [3T, we may think of B as being Haar distributed in ILJS'p(n). Let the first 


quaternionic column of B^ be [gi,..., g„]^, qj G Hr, with 


if- 

Now [/™ = BD^B^ and a simple calculation imply 


{e.XiUre,) = ^ |a,f |/3,f e 


— imOi 


7 = 1 


7 = 1 


(e2,£(C/)"®e2) = E l/S^f e®”®®^+E l«7f e 


— imO-i 




{e^,€iUre2) = -J2 


7 = 1 


«7^7e®-®.+Ed,^, 


7 = 1 


7 = 1 


(e 2 , £(17)"®ei) — - ^ Lx.jyj 
7 = 1 


ai/li,e®”®®" -I- ^ aj/3j 








Therefore the matrix-valued measure (2.7) has the form (3.14) with 

|2 

|2 


IT, = 


a,- 


. *^ 7/^7 1/^7 1 

Using Proposition 2.2 ’c) and Lemma Qb), we obtain the promised joint distri¬ 
bution. □ 
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4. Matrix Models of Unitary Random Matrix Ensembles 


It turns out that if one applies the CMV-fication procedure described in Sec¬ 
tion [ 2 ?^ then the Verblunsky coefficients become statistically independent with the 
distributions that are explicitly computable. For CUE{n) and SO{2n) this was 
done by Killip-Nenciu in [ 22 ]. We restate these results here for completeness and 
develop analogues for the other ensembles. We stress that for the quaternionic en¬ 
sembles CSE and USp, one needs to consider the 2 x 2 block CMV-fication if one 
hopes the same phenomenon to happen. 


Definition 4.1. (a) We say that a real-valued random variable, X, is beta- 

distributed with parameters s,t>0, which we denote by X ^ B{s,t) if 

E{/(W)} = r(s)r(tr^ / ^ -f xy-^dx. 

We write X ^5(0, 0) when X takes values 1 and —1 with probabilities 1/2. 
(This is the t = s — >■ 0 limit of the above.) 

(b) We say that a complex random variable, X, with values in the unit disk D 
is Q{v)-distributed (for v > 1) if 


E{f{X)}=^ I 

where dfz is the Lehesgue measure on D. As a limit case, X ~ 0 ( 1 ) means 
that X is uniformly distributed on the unit circle. 

(c) We say that a 2 x 2 matrix-valued random variable, X, is T(iy)-distributed 
(for iy>5) if 

^ ^ ao -I- aii -02 -I- a^i 

02 -I- 03! Oq — Oil ’ ^ ' 

where, writing a = (oi, 02,03, 04), 


<1 


/(ai')(l — Ix^y'^ ^^^^dxidx2dx3dx4. 

T( 3 ) when X is Haar distributed in SU(2), 


E{/(a)} = / 

J\\s\ 

By extension, we write X 
that is, (4.11 with (01,02,03,04) uniformly distributed over {a; G 
xl -\- xl + = 1 }. 


xj 


Remarks. 1 . Incorporating our stated extensions, B{v(2, v(2) is the distribution of 
the first component of the random vector uniformly distributed over the u-sphere (if 
u > 0). 0(u) is the distribution of a:i-|-ix2, where xi, 2:2 are the first two components 
of the random vector uniformly distributed over the u-sphere (if u > 1 ). Finally, the 
T(u)-distributed 01,02,03,04 appear as the first four components of the random 
vector uniformly distributed over the u-sphere (if u > 3 ). This is shown in [221 
Cor A. 2 ]. 


Proposition 4.2. The CMV-fication of ensembles COE{n) or CUE{n) produces 
the CMV matrix C{aQ,... ,otn-i) with independent Verblunsky coefficients 


Oik ~ Qiffin — 1 — fc) -I- 1), 0 < fc < n — 1, (4.2) 

where /3 = 1 for COE{n) and f = 2 for CUE{n). 

More generally, for any 0 < /3 < 00 , the spectral measure 0 /C(ao,..., a„_i) 
with independent Verblunsky coefficients (4.2) has the law (3.4). 
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Remarks. 1. The ensembles of CMV matrices with (4.2) for any 0 < /3 < oo were 
introduced by Killip-Nenciu in [55], and are called the circular /3-ensembles. One 
way to look at it is that Dyson’s ensembles embed themselves into a family of ran¬ 
dom matrix ensembles that depend on 0 < /3 < oo continuously. This looks natural 
if one thinks of the eigenvalue point process as an interacting particle system. Sim¬ 
ilarly, we will see a /3-extrapolation for SO(n) and 0(n) \ SO(n) ensembles below 
as well. It would be interesting to find a natural /3-extrapolation for USp{n). 

2. The CMV-fication procedure requires cyclicity of ei or, equivalently, that 
\ak\ < 1 for fc < n — 2. This holds with probability 1. 

Proof. The result for CUE{n) was established in [551 Prop 3.3]. 


The result for COE(n) is the combination of our Proposition 3.1 above, [551 
Prop 4.2], and the fact that two CMV matrices with the same spectral measures 
must coincide. 

The last statement of Proposition is in [55) Prop 4.2]. □ 

Proposition 4.3. The 2x2 block CMV-fication ofCSE{n) produces the block CMV 
matrix C(ao/ 2 , • ■ • with independent matrix-valued Verblunsky coefficients 

aol 2 , ■ ■ ■, Oin-iR, where aj are distributed according to ( [4.2] ) with /3 = 4. 

In other words, the 2x2 spectral measures ofCSE{n) and o/C(ao/ 2 , • ■ • > <t„_i/ 2 ) 


with aj given by (4.2) with /3 = 4 coincide. 


Remark. Scalar CMV-fication of CSEin) is not natural. It would produce CMV 
matrix with 2n Verblunsky coefficients with non-independent distributions. 


Proof. By Proposition 3.1 the random 2x2 spectral measure of a CSEfn) ma¬ 
trix is equal to I 2 times the scalar measure (]3.1[) with (3.4) and /3 = 4. This 


implies that the 2x2 Verblunsky coefficients of CSE{n) are equal to I 2 times the 
scalar Verblunsky coefficients of the scalar measure, which are (4.2) by the previous 
proposition. □ 

Proposition 4.4. The CMV-fication of the 0(/V) ensemble produces the CMV 
matrix C{ao,... ,aN-i) with independent Verblunsky coefficients distributed as 

afe 0<k<N-l. (4.3) 

Remark. CMV-fication of SO(n) and 0(n) \ SO(n) gives the same (4.3) except for 
the last coefficient a„_i which has to be taken always —1 for SO(2fc) or 0(2fc-|- 1) \ 
SO(2k -|- 1) or 1 for SO(2fc -|- 1) or 0(2k) \ SO(2fc). 

Proof. This was proved for SO(2n) in [55] Prop 3.5]. The other three cases can be 
proved in the exact same way. Finally note that the total measure of SO(V) and 
0{N) \ SO(/V) with respect to the Haar measure on 0(/V) are equal. This gives 
the 33(0, 0) distribution for otv-i- Q 

The orthogonal groups treated in Proposition |4.4| can be realized as special cases 
of the following more general result: 

Proposition 4.5 (Real orthogonal /3-ensemble). For any 0 < (3 < 00 , a > —1, and 
b > —1, we have: 

(a) (TV = 2n,det = 1) The spectral measure ofC{ao, •. ■ ,a 2 n- 2 , —1) with inde¬ 
pendent Verblunsky coefficients 


Otk 


jd -\- a-\-l, ^ *’ /3 -1-6-1- 1 ), 0 < k < 2n — 2, k even, 

j3 -\- a -\- b -\- 2, ~^ /3) , 0 < k < 2n — 2,k odd 


TZ) { 2 n— 2 —k 

a ( 4 H -r u- -r J., 4 

^ ( 2 n— 3—fc R \ \ h { 2n—l — k 
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is (3.6| on (3.7), (3.8) with the joint distribution 


c„(, 


1^(2 cos 01 ,..., 2 cos0„)|^ IT ~ cos 0^1““'’2 |1 + cos 0^1^“''2 d0i ... d9„ 


3 = 1 


i=i 




(b) (N = 2n, det = — IJ The spectral measure ofC{ao ,..., a 2 n- 2 , 1) with inde¬ 
pendent Verblunsky coefficients 


afc ~ B , 0<fc<2n-2 

is (3.9) on (3.10[), (|3.11 ) with the joint distribution 


(4.4) 


n—1 


D 


. Mi Mi 

2py |A(2cos0i,... ,2cos0„_i)|^ |1 —cos0j| 4 2 |i + cos0j|4 2d9i...d9„ 


1=1 


n —1 


^ ^ IT dp-l---dp,„. 

1 = 1 

(c) (N = 2n + l,det = 1^ The spectral measure of C{ao,..., a 2 n-i,^) with 
independent Verblunsky coefficients 

5(^/3,^/?), 0<fc<2n-l 

is (3.12) on (3.7), ( |3.11 ) with the joint distribution 

M 1. 

1 — cos 9j I 

1=1 


3/S 1 /3 1 

^ |A(2cos0i,..., 2 cos0„)|^ IT ~ cos 0 j|T “2 |l + cos0j| 4 “2 d9i... d9n 




^ M„(0)Tn+l IT 7 */ dfj,i . . . dp,n- 
1 = 1 

(d) (N = 2n + 1, det = — 1^ The spectral measure of C{ao,... ,a 2 n-i,—^) with 

independent Verblunsky coefficients 

0<k<2n-l 


is (3.13) on (3.7); (3.11) with the joint distribution 


E, 


^ 0 1 3/3 1 

|A(2cos0i, ... ,2 cos0„)|^ IT ~ cos0j| 4“2 |l + cos 0 j|T “2 d9i... d9„ 


1=1 


4-1A 4-1 


^ M„{P)dn+l IT 7 */ dfj,i...dp,n. 
1 = 1 
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The normalization constants are: 

^|2n(a+h+l)+/3n(n-l) W r(a + 1 + §j)r(6+ 1 + §j)r(§(j + 1)) ^ 

j=o r(a + & + 2 + |(n — 1 + j))r(|) 

[nm\ 

r(i/3n) ’ 

1 r(fj) 

in(/3) 2(2»-1)|- r(i + fi)’ 

r(i^n) 

1 n!^" r(fj) 

[r(i/3)]"r(i/?) 
r(i/3{n+i)) 

Remarks. 1. If /3 = 2, a = 6 = —1/2 in (a), and /3 = 2 in (b), (c), (d), then one 
obtains matrix models for SO(2n), 0(2n)\SO(2n), SO(2n +1), 0(2n +1)\SO(2n + 
1), respectively. 

2. The ensemble of matrices in (a) was introduced in [52]. It is now called the 
real orthogonal /^-ensemble (see Forrester jS] Sect 2.9]). 

3. This admits a log-gas interpretation. Indeed, the eigenvalues of (a) can be 
viewed as n unit charges , ■ ■ ■, e*®" lying on the half-circle 0 < 9 < tt that also 
interact with their image charges at , ■ ■ ■, and with the fixed charges 
— I -I- at 1 and — | + at 1. As usual, 0 < /3 < oo plays the role of the 
inverse temperature. This is worked out in jSJ Prop 2.9.6]. Similarly, one can check 
that (b) corresponds to the charges 1 — ^ at 1 and —1; (c) to the charge 1 — ^ at 
1 and —^ at —1; (d) to the charge —^ at 1 and 1 — ^ at —1. 


Cn{a,b,l3) = 
KniP) = 
DniP) = 

LniP) = 

Eni/3) = 

Mn{P) = 


Proof, (a) was proved in [551 Prop 5.3]. 

Let us prove (b). First of all note that C is a real orthogonal matrix. Equa¬ 
tion (2.6), Lemma [B.ll i), and Q! 2 n-i = 1 implies detC = —1. Therefore C has 
eigenvalues ,..., ,1,-1. 

Next we note that there is a one-to-one map between 2n x 2n CMV matrices 
of the form C(ao,..., a 2 n- 2 ,1) with all aj G (—1,1) and spectral measures of the 
form (3.9) with ( [3.10 ), (3.11), and distinct 9j’s. Indeed, all Verblunsky coefficients 
are real if and only if the spectral measure is symmetric. And as we just saw, in 
this case a 2 n-i = 1 is equivalent to detC = — 1, which for orthogonal matrices with 
simple spectrum is equivalent to having an eigenvalue at each 1 and —1. 

This map therefore induces a one-to-one mapping between (ag,..., a 2 n- 2 ) G 
(-1, l)2"-i and (9i,.. .,9n-i,Mu ■ ■ ■, 9n) in ([3.10|), (|3.11[) with di < ... < 6»„. So 


we can speak of the corresponding Jacobian. In fact, 
and |3.2[ b), we obtain its value: 


combining Propositions |4.4| 


det 


d(ao, . . . , Q;2n-2) 


d(9 


1 ,... 


: 9n—l 5 Ml 5 ■ • ■ 5 Tn') 


|A(2cos6»i,..., 2 cos6>„_i)|^ IlLi sin^ 


= 2 


n—1 


2n—l—k 


\J MnMn-l-l 11^=0 ^ 
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Using the form (3.9) of the spectral measure, Lemma B.l v) gives 

2 n-2 


na 




j=0 


n—1 

= P ^ |A I 


„ — i9i 


A9ri—1 ^ — 'i'9n — 




i=i 


n—1 n—1 

= II II I sin6»j|^ |A (2cos6*i, ... ,2cos6»„_i)|^ . 

j=i i=i 


Now, starting with the distribution (4.4), and using the change of variables 
along with the last computation produces the promised distribution of the spectral 
measure. The normalization constant Tn(/3) is a well-known integral (see, e.g., [22l 
Lem A.4]). 

Parts (c) and (d) are similar. □ 


Proposition 4.6. The 2x2 block CMV-fication of the I[J5'p(n) ensemble produces 
the block CMV matrix C{aQ,... with independent 2x2 Verblunsky coeffi¬ 

cients Uj distributed as 

otk ~ T(4n — 4A: — 1), 0 < fc < n — 1. 

Remark. Scalar CMV-fication of l]Sp{n) is not natural and does not lead to sta¬ 
tistically independent Verblunsky coefficients. 

Proof. With the help of our Propositi on |2.2K c), the proof for ILJS'p(n) goes along 
the same lines as the proof for l[J(n) in [22[ Prop 3.3]. □ 


5. Matrix Models for Truncations 
In order to deal with the truncations, we will need the following lemma. 

Lemma 5.1. Consider a unitary CMV matrix C(ao,..., ocn) with an G 9D. 

• If n is even then the matrix obtained by removing the first row and column 
is equal to 

( an—ian^ C^n—20^n: ■ ■ • ; S 

for some unitary S. 

• If n is odd then the matrix obtained by removing the first row and column 
is equal to 

C ( an—ian.) C^n—2C^n: ■ • ■ ; OiQan} S 
for some unitary S. 

Proof. Suppose first that n is even. Then C(ao,...,a„) with the first row and 
column removed is equal to CM, where 

C = diag ([-tto], S[q! 2], ■ • •, S[a„_2], [d„]), 

M = diag(S[ai],S[a 3 ],...,S[a„_i]), 


where ^[a^] is defined as in (2.1). 
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Now consider {Q^P 2 M'^pIq){Q^P iC'^P^Q), where 


Q = 

and Pi = diag(a„, 1, a„, 1,...), P 2 
Direct calculations yield 


0 ••• 0 1 

0 ••• 10 

■ . • • 5 

_1 • • • 0 0 _ 

= diag(l,a„,l,a„,...)- 


P 2 M^pIQ = diag(^[-a„_ia„],... ,^[-aian]), 
Q^PiC^pIq = diag([l],S[-Q!„_ 2 a„],... ,S[-d 2 Q;„], [-doan]) 


This shows Q'^P 2 {CM)'^P^Q = P 2 M'^pIQ){Q'^P iC"^P^Q) is precisely equal to 

Cj O-Ti—— 2^715 ■ ■ ■ 5 cf. Definition 

Similarly, if n is odd, we can see that {Q^P^CP 2 Q)(Q^P 2 MP 1 Q) is equal to 
C ( Q-ji—iCym C^n—2C^m ■ ■ ■ 5 Q^O^n)- ^ 


With this lemma and the results of Section in hand, we immediately obtain 
matrix models for the truncations of unitary ensembles. 

Indeed, choose a matrix U from one of the “scalar” ensembles from Section 
{CUE, COE, O, SO, O \ SO), and remove any row and the corresponding column. 
By the invariance under the group action, it does not matter which row and column 
we delete. It will be convenient to delete the first ones. Then the truncated matrix 
V is 

V = PUP\ where P = proj ((ei)"*‘) . 


Recall now (2.5) that U = WCW^ with Wei = W^ei = ei which implies 


V = PWCW^P^ = WPCP^WK 


Therefore V is unitarily equivalent to the CMV matrix form C oi U but with the 
first row and column removed. Thus the results of the next proposition follow from 
the CMV models in Section]^ the previous lemma, and the rotational (or real even) 
symmetry of the Verblunsky coefficients. 

For quaternionic ensembles CSE and USp one needs to modify all of the argu¬ 
ments in the obvious way to accommodate quaternionic rows and columns. 


Proposition 5.2. 

(a) The eigenvalue distribution of a COE{n + 1) or CUE{n + 1) matrix with 
one row and column removed coincides with the eigenvalue distribution of 
C(ao) • ■ •) C(n-i) with independent Verblunsky coefficients 

otk ~ &{P{k -I-1) -|- 1), 0 < fc < n — 1, (5.1) 


where /3 = 1 for COE and f3 = 2 for CUE. 

The eigenvalue distribution of a CSE{n 1) matrix with one (quater¬ 
nionic) row and column removed coincides with the eigenvalue distribution 
of C{aQl 2 ,. ■., < 371 - 1 .^ 2 ) where are independent and distributed according 
to (5.1) with fi = A. 
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(b) The eigenvalue distribution of aO{n+l), orSO(n+l), or 0(n+l)\S0(n+ 
1) matrix with one row and column removed coincides with the distribution 
o/C(ao,... ,a„_i) with independent Verblunsky coefficients 

ak ^ B , 0 < fc < n — 1. 

(c) The eigenvalue distribution of aUSp{n + l) matrix with one (guaternionic) 
row and column removed coincides with the distribution of C{ao,..., a„_i) 
with independent 2x2 matrix Verblunsky coefficients 

T (4fc + 7), 0<fc<n — 1. 


Remarks. 1. From the aesthetic considerations, one might prefer to look at the 
above distributions as B{ ^ i(fc+^ 2 )-i ^^ 0(2(A: + 2) — 1), T(4(/c + 2) — 1) for 

the truncations of the compact groups O, U, and U^p, respectively. Note B has an 
extra factor of ^ due to the historical definition of the Beta function. 

2. These CMV matrices are no longer unitary as the last coefficient is no longer 


unimodular (with probability 1), see the remark after Definition 2.3 
3. These models are efficient for simulation purposes. 


6. Eigenvalue Distribution of the Truncations 

6.1. Truncations of circular /?-ensembles. 

Theorem 6.1. Choose Verblunsky coefficients ak to be independent and 0(/3(fc 
1) + 1)-distributed (0 < k < n — 1), that is, distributed in D with to the density 




Then the eigenvalues of the CMV matrix C{aQ,... ,an-i) are dis¬ 

tributed in D” with the joint law 


j,fc=l j<k 


( 6 . 1 ) 


Remarks. 1. For any 0 < /3 < oo these ensembles are unitarily equivalent to the 
circular /3-ensembles with the first row and column removed. 

2. Combining this with Proposition |5.2[ we obtain Theorem |1.2| 

3. As we will see in Section]^ the law (6.1) coincides with the Gibbs measure 
for a system of n charges in the disk D where D and C \ D are filled with dielectrics 
of suitably chosen susceptibility. 

The proof is given in the end of this subsection after we establish a series of 
lemmas. 

Given a sequence of Verblunsky coefficients uq, ... one may then form 

the associated sequence of orthogonal polynomials $ 0 )--■)'!>„; see Appendix [B| 
Recall from (2.6) that the eigenvalues of C(ao, • ■ • ><Tn-i) coincide with the zeros of 
$„(z), which we shall denote by zi,... ,Zn. Note that the mapping from the zeros 
{zi,... ,Zn) € D” to the Verblunsky coefficients (ao,... , Un-i) G D” is an n-fold 
branched covering. This follows from 1301 Thm 1.7.5]. To compute the Jacobian 
of this change of variables, it will be useful to introduce an intermediate set of 

(k) 

variables. Let n) be the coefficients of the fc-th orthogonal polynomial <l>j,: 
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The Jacobian of the mapping from roots of a polynomial to its coefficients is 
well known; we reproduce it here as Lemma Qb). Viewing complex numbers as 
ordered pairs of real numbers, we have 


9(4 


, Kn 


0(^Z\ , . . . , Zfi') 

Lemma 6.2. The Jacobian Ji 






is equal to (-l)”n4o(l“ 1“*^) 


(6.3) 


\2\k 


Proof of Lemma 6.2. From (B.21 we can immediately see that 


(n) ('^—1) - -(^—1) • 1 

4 j = 1,... ,n, 


(6.4) 


where we assign Kn" = 0 and Kq" = 1. Therefore applying Lemma 
obtain 


D.l 


we 




d{. 


.(—1) 




,-i) 


det 


r 1 

0 

0 

Cln— 1 

0 

0 ■ 

• O 

1 

^n—1 

0 

0 

0 

0 

Hn— 1 

1 

0 

0 

0 

^n— 1 

0 

0 

1 

0 

0 

0 

(n-1) 

-<-l •• 

0 

(ra-1) 

-«1 

0 

-1 

-(n-1) 

L-<-i 

0 

-(n—1) 

• -«:i 

0 

-1 

0 . 


. (6.5) 


Expanding this determinant along the last two columns and then rearranging the 
remaining rows and columns we obtain a block diagonal matrix with n — 1 blocks 


1 

1 


^n — 1 

1 


on the diagonal. Therefore (6.5) simplifies to —(1 — |q;„_iP)"“^ (minus sign comes 
from expanding along the last two columns). 

The lemma now follows by induction. □ 


Lemma 6.3. The Jacobian 


d(ao. 


.-i) 


a(2i,...,z„) 


is equal to 






Proof. This is immediate from (6.3) and Lemma 6.2 


□ 


Now we can prove Theorem |6.1[ 


Proof of Theorem \6.1\ The joint distribution of ag,..., an-i is 

n— 1 

St n (1 - 

fc=0 
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Using Lemmas 6.3 and |B.l[ vi) we see that this induces the following measure on 


/n—1 


n! ( 27 r)^ 




\k^0 






fe=0 


= j£^\A{zi, . . . , Zn)\‘^ ]J( 1 -Zj 0 fe )2 (fzi . . .(fzn- 
j,k=l 


Note that we removed the ordering of Zj^s which resulted in the ^ factor. □ 


6.2. Truncation of the orthogonal /3-ensemble. Truncations of the orthogonal 
group, and more generally of the orthogonal /3-ensemble, have real eigenvalues with 
non-zero probability. Correspondingly, we decompose the set of possible systems of 
eigenvalues as 

D := [J Dl^m, 

L+2M=n,L>0,M>0 

where 


Dl,m ■= {(xi,..., xl,xl+i ± iyL+i, ■ • ■, xl+m ± Wl+m) € : 

-1 < Xi < X2 < ■ ■ ■ < XL < ^, -1 < xl+1 < ... < Xl+m < 1, 

j/L+i > 0,... ,yL+M > 0} . (6.6) 

Assume we are given a function / : C" —)■ C" that is non-negative on D and in¬ 
variant under the permutations of its n arguments. Let us define f{zi ,..., Zn)\dziA 
dz 2 A ... A dzn I to be the distribution on D C D" determined by 

/ f{zi,... ,Zn)\dzi AdZ 2 A ■ ■ ■ AdZn\ 

JB’i 

._ ^ 2^ f{xi,...,XL,XL+lHyL+l,---,XL+MayL+M) 

L+2M=n dDL,M 

X dxi... dxL+ndyL+i ■ ■ ■ dyL+M- 

Informally one should think of each of the summands on the right hand side as the 
distribution conditional on the event that the matrix has L real eigenvalues and M 
complex-conjugate pairs with L + 2M = n. Each of these distributions is obtained 
by plugging into the left-hand side dzk = dxk + idyk, taking the wedge product 
while keeping track which of the eigenvalues are real or complex, and then taking 
the absolute value of the constant. Thus for each real eigenvalue Zk-, dzk becomes 
just the real Lebesgue measure dxk on the real line, and for a complex-conjugate 
pair Zk, Zk, we get \{dxk + idyk) A {dxk — idyk)\ = ‘^dxkdyk- This corresponds to 
the factor 2^ on the right-hand side of the definition. 

More formal and careful discussion of these measures can be found, e.g., in the 
Borodin-Sinclair paper [5J Sects 2-3]. 
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Theorem 6.4. Given /3 > 0, a,b > —1, let the coefficients ak, 0 < k < n — 1, be 
distributed as follows: 


ak 


B 


B + a + 1, ^ + 6 + , if k is even, 

(^^^+a + b + 2,^^y ifkis 


odd. 


Then eigenvalues of the CMV matrix C{ao ,..., On-i) are distributed in . 
joint law 


(6.7) 


with the 


jr Ylil-ZjZk)' ^Y[il-ZjT+^ 

j,k=l j = l 

where 

L^J 

p _ 2"(“+f>+i)+f 1) 


^ \Zk-Zj\ \dz1Adz2A.. .AdZn\, 
j<k 


r(a + 1 + ^_))r (6 + 1 + ^j) 
r(a + b + 2 + I (LfJ + j)) 


LtJ 


nr(§j)- (6.8) 


^ ^ , 2 VL2J ' -'ll j = l 

Remarks. 1. From our previous discussion it is evident that if n is odd, then 
such a random matrix is unitarily equivalent to a (n + 1) x (n + 1) matrix from the 
orthogonal /^-ensemble in Proposition |4.5[ a) with the first row and column removed 
(for any 0 < /3 < 00 and a,b > —1). If a = & = —1 + /3/4, then it is also unitarily 
equivalent to a (n + 1) x (n + 1) matrix from the ensemble in Proposition 4.5 b) 
with the first row and column removed (for any 0 < /3 < 00). If n is even and 
a = b — —1 + /3/4, then it is unitarily equivalent to a (n + 1) x (n + 1) matrix 
from the ensembles in Proposition |4.5[ c) and (d) with the first rows and columns 
removed. 

2. When ft = 2,a = b = —1.12,, such a random matrix is unitarily equivalent to 
a matrix from 0(n + 1), or SO(n + 1), or 0(n + 1) \ SO(n + 1) with any row and 
column removed, cf. Proposition |5.2K b) . 

We follow the same strategy as for the circular ensembles. As before, let (6.2). 
In the next lemma we compute the Jacobian determinant of the transformation 
from Zj’s to on each of the sets Dl,m- 


Lemma 6.5. Suppose Zi,..., Zn are in Dl,m, L+ 2M = n, see (6.6\. On this set 
we have that the mapping ( 21 ,..., z„) 1 —>■ (ki"\ ..., k^T^) is injective, and 


J 




.. . =2"A(zi,...,2„). 

0{xi, . . . , XL, Xi+I, l/L+l, a;L+2, J/L+2, ■ • • ,XL+M,yL+M) 

Remark. In other words, the above lemma says that the following change of vari¬ 
ables is legitimate 

/(...) dn^f'K .. dnlf-'* =/(...) |A(zi,...,z„)| \dzi A dz 2 A ... A dzn\. 


Proof. Note that the above change of variables can be viewed as a composition of 
two changes of variables: 

(«:^”\ ..., !-)► (xi,..., XL, XL+i ± ipL+i, • ■ ■, xl+m ± Wl+m) 

!->• (xi,. . . , XL, xl+ 1, J/L-l-1, a:L+2, yL+2, • ■ • , Xl+m, yL+M) 

But the (complex) Jacobian determinant for the first map is A(zi,..., Zn) by 
Lemma [D.2[ a), and the (complex) Jacobian for the second map is easily seen to be 
2^ in the absolute value. □ 
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Lemma 6.6. The Jacobian determinant of the 


, K^n’) 

d{ao,..., a„_i) 


—>■ M" change of variables 

=(-ir n n 


0<fc<n-l 
k even 


0<fc<n—1 
k odd 


Proof. Using the recurrence (6.4) it is easy to see that 


J 




= det 


1 

0 

0 

1 

(n-1) 
n-1 


C^n—1 
0 


— K. 


— K. 


>-l) 

n-2 


0 

O^n— 1 

1 

0 

_ ("-!) 
A^n 


O^n—1 

0 


0 

1 

d"-l) 


0 ■ 
0 

0 

0 

-1 


(6.9) 


If n is even then the (^, ^)-entry of the latter matrix is 1 — an-i- We stress that 
this is an n X n matrix, in contrast to (6.5) which was 2n x 2n. Expansion along 
the last column allows to compute (6.9), and then induction finishes the proof. □ 


Now we are ready to prove Theorem 6.4 


Proof of Theorem \6.4\ Starting with the joint probability density (6.7) of a^’s and 
performing two changes of variables from the previous two lemmas, we obtain 


7^ n (1 “ “fc) 


_ -l-n _— 

4 0 . 


Hl + ak)^ 


A 9 


0<k<n-l 
k even 


n (fc-l)/3 I I L I -I fc + 1 , , (fc + l)/3 .. fc-1 

(1-afc) 4 +“+6+1 2 (1+Q,fe) 4 1 2 


0<A:<n— 1 
k odd 


where the constant 

n 


1 

Pn 


X |A(zi,... ,z„)| \dzi A (iz2 A ... A dz„|, 

r(ffc + a + 6 + 2 ) 


0<fc<n—1 
k even 


2T+“+6+ir(f A: + a + l)r(f fc + 6+1) 


n 


V(^^k a b 2) 


0<k<n 
k odd 


2^+"+^+ir(f (A: - 1) + a + 6 + 2)r(f (A: + 1)) 


can be easily seen to be equivalent to (6.8). 

After some simplification, the distribution becomes 

n—1 n—1 n—1 

7 + n (1 - n (1 - n 


^ A^ 9 






k^O 


X |A(2:i, ..., Zn)\ \dzi A dz2 A . .. A dzn\. 

Now, an application of Lemma [B.l|[i) and (vi) finishes the proof. □ 
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7. Matrix Models and Resonance Distribution for Systems with 

Non-ideal Coupling 


As argued by Fyodorov-Sommers in [TS] (see also [HIIS]), the eigenvalues of 
the truncations represent resonances for the open systems for the special case of 
the ideal coupling. They proposed a more general model (for CUE) which for the 
1 open channel case takes form 

y = C/diag(Vl-7;,l,l,...,l). (7.1) 

Here 17 is a unitary random matrix, and 0 < Tq < 1 is the so-called transmission 
coefficient, which may be random, or deterministic. The case of Tq = 0 corresponds 
to the closed system, and Ta = 1 to the open system with the ideal coupling. Clearly 
the eigenvalues corresponding to Tq = 0 and Tq = 1 (ignoring the zero eigenvalue) 
coincide with the eigenvalues of the unperturbed unitary random matrix and of its 
truncation, respectively. 

In fact, instead of the transmission coefficient it will be more convenient for us 
to work with the reflection coefficient Ra := ^/l — Ta. 

Our method of reducing to the CMV form applies here as well. Indeed, let U be 
n X n and its CMV form he U = WC{aQ, ..., a„_i)IF'^ with Wci = W^Ci = ei. 
Then 

V = WC{ao, a„-i) diag {Ra, 1,1, ■ • ■, 1) 

Using the idea in Lemma |5.1[ one then obtains that V is unitarily equivalent to 

C ( CXn—2^n—l, Oin—3^n—l, ■ • ■ , 17 Eae^n—l') 

if n is odd and to 

C ( 0^,a—2^n—l, ^n—3*^n—1; ■ • ■ , 1: 

if n is even. 

Combining this with the CMV models we developed earlier, we immediately 
obtain the following proposition. 


Proposition 7.1. Let 0 < Ra < \ he independent of the random matrix U. 

(a) The eigenvalue distribution of (7.11 with U G COE{n),CUE{n),CSE{n) 
coincides with the eigenvalue distribution of C {ao,... ,an- 2 , RaE‘^) with 
independent coefficients 

Ofc - 0(/3(fc-t I)-h I), 0<fc<n-2, e*‘^G0(I), (7.2) 


where /3 = I, 2,4 for COE, CUE, CSE, respectively. 
(b) The eigenvalue distribution of (7.1) with U G 0(n), 


or 


n), 


or ' 


(n) \ 


(n) coincides with the distribution ofC{ao,..., a„_ 2 , RaCr) with indepen¬ 
dent coefficients 

0 <fc<n-2, 


Tfc 


and a G B{0, 0) for 0(n); a = —I for SO(2fc) or 0{2k -|- I) \ SO(2fc -|- I); 
and a = 1 for SO(2fc -|- I) or 0(2fc) \ SO(2A:). 

(c) The eigenvalue distribution of a USp{n) matrix with one row and column 
removed coincides with the distribution of C{ao,..., an- 2 , RaQ) with inde¬ 
pendent 2x2 matrix coefficients 


ak--T {Ak + 7), 0<k<n-2, g G T(3). 
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Remarks. 1. As usual, CSE has double multiplicity at each of the eigenvalues 
of ( [A2| . 

2. When = Q these CMV models coincide with the models in Section 
(with n instead of n + 1 that we had earlier), and when Ra = 1 they become the 
(reversed) models from Section (the reversed models however do not have the 
spectral measures preserved but only the eigenvalues, so are less natural). 

For (a) and (b) we can compute the eigenvalue distributions. 

Proposition 7.2. (a) Let 0 < (3 < oo and 


<Xk 


0(/3(fc+l) + l), 0<fc<n-2, e*‘^e0(l). 


Assume the reflection coefficient Ra is independent of ao,... ,an-2,4> with a dis¬ 
tribution on (0,1) of the form 


F{Ra)fn{l-Rl)^'^-^RadRa 


(7.3) 

where F is any function that makes ( |7.3| ) a probability distribution. Then the eigen¬ 
values zi,... ,Zn ofC{ao ,..., ctn- 2 , RaF'‘^) are distributed in D” with the joint dis¬ 
tribution 


j,k=l 


j<k 


F 


ln^ 

i=i 


d^Zi . . . dfzn. 


(7.4) 


(b) Let 0 < /3 < oo and 
ctk^ B 


/ /3(fc+i) 

' 4 ’ 4 ' 




f 5 


0 < k < n — 2, B{0, 0). 


Assume the reflection coefficient Ra is independent of ao,... ,an- 2 ,o’ with a dis¬ 
tribution on (0,1) of the form 


FtR I 4 "+2 ) ("I 


where G is any function that makes (7.51 a probability distribution. Then the eigen¬ 
values Zi, ... ,Zn of C{aQ, 
tribution 


>cin- 2 , RaO') are distributed in 


^dRa, (7.5) 

ligen- 

with the joint dis- 


P„ ZjZk) ‘ 

j,k=l 


n 

i=i 


1 - Z] j<k 


n \^k-zj\Gi I 


i=i 


\dz\ A dz 2 A ... A dzn\^ 
(7.6) 


where P„ is (6.81. 


Remarks. 1. When /3 = 1,2,4 in (a), (7.4) is the eigenvalue distribution of (7.1) 
with U e GOF{n),GUF{n),GSE{n), respectively. Also, for any 0 < < oo, (7.4) 

is the eigenvalue distribution of (7.1) with U a matrix from the circular /3-ensemble. 

2. When /3 = 2 in (b), ( |7.6[) is the eigenvalue distribution of ( |7.1[ ) with U G 0(n). 
The analogue of Theorem |6.4| for general a,b > —1 can also be stated. We chose 
a = b = —1 -\- j3/4 to simplify the presentation. 

3. When F = 1 in (a) (or G = 1 in (b)), we obtain that a system under the 
non-ideal coupling with such a reflection coefficient Ra behaves as a system one 
size larger with ideal coupling. One should also recognize that (7.3) with F = 1 
and (7.5) with G = 1 coincide with the distribution of |an-i| in Theorem 6.1 and 
Theorem |6.4[ respectively. Ultimately, this is a manifestation of the statistical 
independence of the Verblunsky coefficients. 
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Proof. The result is obvious from the arguments in Section T he on ly additional 


fact we need to use is that Ra = |q!„_i| = 


m= 




see Lemma 


B.l 


0 - 


□ 


8. Log-gas Interpretation 


In this section, we present an interpretation of the measure ( 6.11 as the Gibbs 


measure for a configuration of charges in a certain geometry of dielectrics. What 
follows was very much inspired by [U Sect 15.8-15.9]. 

Consider a system of n particles of equal charge -1-1 that are confined to lie in 
the unit disk D of the complex plane. The unit disk D is filled with a homogeneous 
(linear) dielectric of permittivity ei, while the rest of space C \ D is filled with 
a homogeneous (linear) dielectric of permittivity 62 - As usual, one may regard 
two-dimensional point charges as representing parallel lines of charge in a three- 
dimensional setting. 

Recall (cf. §7]) that the electric potential due to a unit point charge at zq in 
the presence of dielectric media is the solution to — Vz • [e(z)VzV(z)] = 6 {z — zq). 
For the configuration described above, and Xq S D the solution is 


V(z\zo) = 




logjl - zzo\ 


2 G D, 
z G C \ D. 


From this, one finds (cf. Problems 3 and 4 in [25l §7]) that the polarization of the 
dielectrics exerts a force on the charge at zq G D of the form 

F = which has potential IT(zo) = - j^^f^logil - jzoP]. 

This force pushes the charge toward the origin when ei > €2 > 0 and toward the 
boundary in the case e 2 > ei > 0 . 

The above calculations reveal that the total energy of a system of n charges 
located at the points C D is then 

H = J2wiz,)+^^Y.V{zk\zj), 

3 

and consequently, the associated Gibbs measure at temperature T has density 


^-H/iksT) 




( 8 . 1 ) 


j,k=l 


3<k 


with exponents 

a = and 7 = 75 —■ 

£1+62 ' lireiksT 

Note that any combination of 7 > 0 and a G K is within the realm of physical 
possibility. In the traditional setting (as elaborated in [5S], for example) the ratio 
£ 2/61 can lie anywhere in the interval [ 0 , 00 ], with the case e 2 /ei = 00 , for example, 
being realizable by using a conductor in the exterior region C \ D and vacuum in 
D. On the other hand, the physical consistency of materials with negative (or zero) 
permittivity has been known for some time. In recent years, such materials have 
even been manufactured; see, for example, for an introduction to the literature 
on this. 


Gomparing (6.11 and (8.1), we see that the distribution of eigenvalues for trun¬ 


cated Dyson models coincides with that of a gas of charged particles in the presence 
of dielectrics; specifically, one chooses 7 = 2 and a = f — 1 . 
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9. Symmetric CMV Matrices 


As we mentioned in Subsection |2.2[ in the case of a system with time-reversal 
invariance {COE and CSE) reducing its scattering matrix U to the CMV form is 
not natural since this breaks the initial symmetry. In this section we introduce a 
symmetric analogue of the CMV form which solves this issue. We will call these 


matrices (9.11 “symmetric CMV”. 


In what follows let ReT = f (T -I- Tl) and ImT = ^{T — T^). 

Proposition 9.1. Let U be a unitary operator on with cyclic vector ei and 

spectral measure (2.3). Then applying Gram-Schmidt orthonormalization procedure 


to ei, [Im{7]ei, [Ret/]ei, [Im(t/^)]ei, [Re([/^)]ei,... produces a basis {s/c} in which 
U has the form 


S (oq, oi,...) :— Af, 


where L = lE is as in (2.2), and Af is the unitary matrix 


A/' = diag([l],4'i,«'3,...), 


= 


V^2(l - Reofe) 


i{l - Ofc) -ipk 
Pk 1 — ttfc 


(9.1) 

(9.2) 

(9.3) 


The matrix S is 7-diagonal and symmetric S"’" = S. In fact, U = RSR^ with 
Rf R = RRf = I, Rei = R^ei = ei. Moreover, if C/^ = U, then R is orthogonal. 


Proof. By the spectral theorem (cf. (2.4|), we may conflate U with the operator 
f{z) I—> zf{z) in L^{dp) and ei with the constant function 1. In this way, we identify 
the basis {sfej^g with the functions formed by orthonormalizing the sequence 
{I,sin0,cos0,sin20,cos20,...} in L^{p). Here and below we write z = e*®. 

Following the notation from m Section 4.2], let us write {x/c}^o 
for the orthonormal bases formed by orthonormalizing {1, z, z~^, z^, z“^,...} and 
{1, z“^, z, z“^, z^,...}, respectively, in if {dp). Recall that the traditional CMV 
matrix C = CAA represents U in the {xk} basis and that AA represents the change 
of basis from {xk} to {xk}; specifically, AAjk = {xj,Xk)- 

It is elementary to verify that Af^N = AA and consequently, 

AfCAf-^ = NCAAN-^ = NCN^. 

Thus, the key assertion to be proved is that N represents the change of basis from 
{Xk} to {sk}. To this end, let us first observe that 

X2n-l{z) = '^4>2n-l{z) = pip 2 ..^p 2 „_i [P2n-1^" H-1" Oz "] , 

1 


X2u{z) = Z ”((';„(z) = 


— [— 0 : 2 n-lZ^ Z "] . 


, z^ "}; we will continue 


Here ellipsis indicates terms in the linear span of {z"”^, 
this convention for the remainder of the proof. 

Proceeding directly from the definition of the Gram-Schmidt process, we have 


S2n-l{z) = C- 


2 % 


0 


2;”’+2 


for some positive constant c. Correspondingly, there is a c' > 0 so that 


S2n-\{z) = 


-y/2(l — Rea2n-i) 


[iP2n-lX2n{z) - i{l - Q; 2 „-l)X 2 n-l] + 
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In fact, (/ = 1 and the ellipsis terms vanish. To see this, we observe that S 2 n-i, 
X 2 m and X 2 n-i are all unit vectors and are all perpendicular to elements in the 
span of ..., This confirms all odd-indexed columns of . 

To verify the even-index columns (after the zeroth, which is trivial) we first 
observe that there are some 7 > 0 and 6 G C so that 

S 2 uiz) = + • • • + 7 ^^^- 

Noting also that 




— 2(1 — Rea2n-i) 


2 +2 


2 (Im a 2 Ti-i) - 


we deduce that 


S2niz) = 


\/2(l — Rea2n-i) 


[(1 — a2n-l)X2n(^) +P2n-lX2n-l(^)] +6' S 2 n-l{z) + ’ 


for some 7 ' > 0 and <5' € C. Again the ellipsis terms vanish because they are 
perpendicular everything else that appears in this formula. To determine 7 ' and 6', 
we first exploit that S 2 n is perpendicular to S 2 n-i. Using the expansion of S 2 n-i in 
terms of {xfc} proved above, this shows that 6' = 0. Lastly, the fact that S 2 n is a 
unit vector confirms 7 ' = 1. This justifies the even-indexed columns of J\f~^ = 
and so completes the proof of (9.1). 


The symmetry and 7-diagonal properties of S are obvious from the construction. 
Finally, suppose t/^ = U. Note that columns of R are the vectors obtained from 
the orthonormalization of ei, [ImUJci, [ReC/]ei, [Im(C/^)]ei, [Re(C/^)]ei,... But W = 
U (entrywise complex conjugation), which means that each of Re)!/"*) and Im(C/"*) 
are real matrices for any m. This implies that R is also real. □ 


Remarks. 1. Similar arguments show that, unless n = 2k and det U = —1 (which is 
equivalent to a„_i = 1 ), an n x n unitary matrix U with ei cyclic can be reduced 
to the “symmetric CMV” form 

S{ao,.. .,an-i) = NCN'^, 

where U = diag ([1], 4'i, ..., ^'„_ 2 ), C = diag (Sq, 52 , S 4 ,..., S„_i) if n = 

2k + I and Af = diag ([ 1 ], ^'i, ^' 3 ,..., T„_i), £ = diag (Sq, S 2 , 54 , • ■ •, S„_ 2 ) if 
n = 2k. As before, 5„_i = while '^n-i = [yon-i], where the square root 

corresponds to the argument Argjg 2 ,r) with the branch cut along ]R_|_ (this is just 
the (1, l)-entry of (9.3) using |a„_i| = 1). 

2. Similarly, unless n = 2k and det C/ = 1, then one can instead work out 
an alternative symmetric CMV form S form obtained by orthonormalizing the 
sequence {ei, [Ret/]ei, [Imt7]ei, [Re(C/^)]ei, [Im([/^)]ei,...}. In this case we obtain 
the matrix 

5 ( q ! o , cti,...) := AfCAf^, 


where £ is as in ( 2 . 2 ), and 


A/' = diag(^[l],Ti,«' 3 ,...j , 

1 + dfe 


'Sk = 


y^2(l -|- Re ak) 


ipk 


Pk 

-i{l + ak) 
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Note that det U = —1 for COE{n) is a probability 0 event. Therefore with prob¬ 
ability 1 any matrix from COE(n) can be orthogonally reduced to the symmetric 
CMV form, while fixing e\ and preserving the spectral measure. Moreover, matri¬ 
ces with different spectral measures lead to different sets of Verblunsky coefficients, 
and therefore different symmetric CMV forms. 

The infinite unitary CMV matrix with = 0 is the bilateral shift (after re¬ 
ordering basis elements). The symmetric CMV form of this operator is 


C2 

0 


0 

i 

^/2 
0 5 

0 


1 

0 


0 0 


0 0 


1 
2 

2 

i 0 
i 0 

0 -2 
0 -I 


0 

0 

0 

1 

2 

0 

0 


Similar arguments lead to the analogue of Proposition 9.1 for 
corresponding (complex quaternionic) symmetric CMV form is just 
as above), and the reducing matrix R is in US'p(n). 


CSE(n). The 
S ^ I 2 (with S 


Appendix A. Quaternions 

Let us carefully define the algebras of real quaternions Hr and complex quater¬ 
nions He, as well as discuss quaternionic matrices. 


A.l. Real and complex quaternions. In modern language, the real quaternions 
were introduced by Hamilton as the unital associative algebra Hr over M generated 
by 1, i,J, k with (non-commutative) multiplication defined by 

= ijk = -1. 

We identify the underlying vector space with by using 1, i, k,J as the standard 
basis. Note the deliberate reordering! This makes for a prettier regular represen¬ 
tation q i-> 2{q): let q = a + b\ + ck + d] £ Hr, then 


q{w -I- xi -|- yk -1- z\) = {w' + x'\ -\- y'k + z']) 

if and only if 


'w'' 


w 


a 

-h 

—c 

-d~ 

x' 

y' 

= 2{q) 

X 

y 

with 2,{q) := 

h 

c 

a 

-d 

d 

a 

—c 

b 

z' 


z 


d 

c 

-b 

a 


(A.l) 


The 2x2 blocks appearing in (A.l) correspond to the standard matrix rep¬ 
resentations of complex numbers. This leads to the usual representation of real 
quaternions as 2 x 2 complex matrices: 


a -I- 6i -I- ck -I- dj !->• 


a-\- ib 
c + id 


—c + id 
a — ib 


<L{a -I- &i -I- ck -I- dj). 


(A.2) 


where a, b,c,d G K. 
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In particular, this shows that the group of unimodular real quaternions (i.e., 
q e Hr with = 1) is the special unitary group SU(2). 


The analogue of (A.l) for multiplication on the right is as follows: 

—b —c 


Dl(a + 6 i + ck + dj) := 


a 

b a 
c d 
d —c 


which is to say 


(■u; + a:i + 2 /k + zj)q = {w' + x'\ + y'k + z'j) 


-d 

a 

b 


-d 

c 

-b 



'w'' 


W 


x' 

n,' 

= in(g) 

X 


y 

y 


z' 


z 


Note that q i-A 9\{q) is not a representation (i.e. homomorphism) of Hr, but 
rather of the reversed algebra. Needless to say, all matrices ^{q), q S Hr commute 
with all £(p), p S Hr. Indeed, the general theory of matrix algebras (cf. [Ml 
Theorem 3.4.A]) guarantees that every matrix that commutes with all £(p), p € Hr, 
must be of the form fH(q) for some q G Hr. 


The complex quaternions He are defined in the same exact way as Hr but over 
C instead of M. For q = a + bi + ck + dj G He (now a, b,c,d G C), we will be using 
the representation €{q) defined by same expression (A.2). 


Note that £(Hr) consists of 2 x 2 complex matrices with an = 0:22 

and ai 2 = —0:21, and £(Hc) consists of all 2 x 2 complex matrices. 

For a (real or complex) quaternion q = a + bi + ck + dj let us define its conjugate 

by 

g = a — 6i — ck — dj, 

and its Hermitian conjugate by 

= d — &i — ck — dj. 

Of course g = g'^ for g G Hr. 


A.2. Matrices of quaternions. For the rest of the section, let Q be an n x n 
matrix of (real or complex) quaternions Q = [qijj^j^i- 

We define £(<5) as the matrix formed by replacing each quaternion entry by the 
corresponding matrix representation 


e:(Q) = [£(?,,)]■,=!. 


Let us define 
It is easy to see that 




the usual Hermitian conjugation of complex matrices. 

For any 2n x 2n complex matrix M we define its (time reversal) dual matrix by 

:= 


where 



-I 
0 ’ 


(A.3) 
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the 2n x 2n block diagonal matrix with n copies 

For a quaternionic matrix Q let us define = 
that if Q is complex quaternionic then 


0 

1 


-1 

0 


along the diagonal. 


£ ^(£((5)^). It is easy to verify 


= fed-,=i. 


If Q is real quaternionic matrix, then coincide, which implies 

Z^€{QfZ = €{Q)^. (A.4) 

We say that a (complex or real) quaternionic matrix is self-dual if = Q. 

We say that a (complex or real) quaternionic matrix Q is unitary if £(<5) is 
unitary. 

Finally, by the eigenvalues of a quaternionic matrix Q we mean the eigenvalues 
of the complex matrix £((5). 


Appendix B. Basics of Orthogonal Polynomials on the Unit Circle 


We will introduce some basics of the theory of orthogonal polynomials on the 
unit circle. For more details we refer the reader, e.g., to mM- 

Suppose we have a positive probability measure /i on the unit circle 9D having 
at least n points in its support. For the purposes of this paper, the typical example 
one should keep in mind is ^ being the spectral measure of an n x n unitary matrix 
U with distinct eigenvalues (see Section]^. 

Then applying Gram-Schmidt procedure to the monomials with respect 

to the inner product {f{z),g{z)) := f{e^^)g{e^^)dfj,{6) in L^(d/i), one obtains a 

sequence of monic polynomials orthogonal with respect to g 

($j(z),d>fe(z)) = 0, (B.l) 


The famous result of Szego is that these orthogonal polynomials satisfy the 
recurrence relation 

d>fe+i(z) = z$fc(z) - dfcd>fc(z) (B.2) 

for some sequence of complex coefhdents {ajljUg from D (referred to as Verblunsky 
coefhcients), where for any polynomial P{z) of degree k we define 


P*{z) = z'^P(z-i) 

(i.e., if P{z) = z'^ + -|- ... -I- Kk-iz + Kk, then P*{z) = RkZ^ + Kk-iz’^~^ + 

KlZ 1). 

It should be noted that if g, has exactly n points in its support, then {z^}^Zq form 
a basis, and therefore z" is a linear combination of {z^}^Iq in L^{dg). However 
Gram-Schmidt procedure can still be applied to produce a polynomial d)„(z) (which 
is of norm 0 in L^{dg), of course), and (B.2) holds with a unimodular coefficient 

CXji— 1 ■ 

In fact, there is a one-to-one correspondence between all probability measures 
supported on n distinct points 


n 

(B.3) 

i=i 

and sequences of Verblunsky coefhcients {ctk} kZo with afcGDforO<fc<n — 2 
and a„_i £ clD. So ao,...,Q;„_i can be viewed as a change of variables from 
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, e*®", ^ 1 ,..., (subject to Mi = !)• This was initially explored 

in [22] (see also Forrester-Rains HD), and is developed further in our Section |4) 
For each 0 < fc < n — 1, denote the orthonormal polynomials 

4^k{z) = ■jjj^$fc(z), 

where the norm is in L^{dii). We summarize some of the properties of orthogonal 
polynomials that we will need in the next lemma. 

Lemma B.l. (i) Let 1 < k < n and (z) = Ilj=iiz- Zj). Then 

k 

i=i 

If all aj’s are real, then 

k k — 1 

®fc(l) = 11(1 - Zj) = 11(1 “ aj); 

j=l j=^0 

n k — 1 

^ fc (- l ) = 11(1 + = 11(1 + {- lyaj ). 

j=l j=0 

(ii) For any Q<k<n, = 11^=0 (1 - Wi?)- 

(hi) For any Q < k < n — 1, the Christoffel-Darboux formula holds: 


j, ( ^X77^ ^kiz)4>k{C) ~ 4’k{z)4>k{C} 

y^yiizmiO = --■ 


j=o 

(iv) The Cauchy formula holds: 

1 


det 


1 - ZiZ. 


n 


.j<s 


LZs - 2 ,- 




ZjZs) 


(v) If fi is of the form (B.3), then 

n—2 


n (1 -=IA ,..., e*^") I n 

j=0 j-l 


(vi) For any 0 < k < n, if ^k{z) = ni=i('^ ~ Zj) then 


n (1 -1«. 

j=0 


= n (1 “ 

j,S=l 


(B.4) 


Proof, (i) is immediate from ( B.2| ). (ii) and (iii) are standard facts, see [301131]. A 
proof of (iv) can be found in, e.g., [52J Lemma 7.6.A], (v) is proved in [221 Lem 4.1]. 
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Let US now prove (vi). An alternate proof of this in [30l Thm 2.1.4]. By elemen¬ 
tary row operations, 




= det [$ 8 - 1 ( 2 :, 


■J)h<j,s<k 


^j-l{Zs) 




/k-l 


KS^O 






Note that 0 ^( 2 ) = — zzj) and (j)k{zs) = 0 for any s. Therefore the 

right-hand side of the last equation can be written as 

= (nV«lP) det 


\s^0 

^k-1 


(t>k{ZjWki^s) 


1 - ZjZs 


= n n \^kizs)? det 


1 


\s^0 


1 


Now using parts (ii), and (iv), we obtain 


|A(zi,...,Zfc)p Ul=i\^-z,zs\^ 


which gives us (B.4| after cancelling by |A| 


□ 


Appendix C. Some Distributions on Spheres 


For the reader’s convenience we collect here some of the distributions needed in 
the body of the text. 

Lemma C.l. (a) Let a random vector x € K" be chosen according to the normal¬ 
ized uniform distribution the unit sphere {x S M” : ||a;|| = 1}. Then 

(^ 1 ,..., fin) — ■ ■ ■, Xn) 

are jointly distributed on the simplex Tj = ^ 0,1 < j < n, according to 

the probability distribution 

n 

WWW n . . . dfin-l. 

i=i 

(b) Let a random vector x G be chosen according to the normalized uniform 

distribution the unit sphere {a; G : ||x|| = 1}. Then 

{fll, • ■ . , fLn) {Xi X 2 , • • • ) X2n—1 T ^2n) 

are jointly distributed on the simplex Lj — 1) Lj ^ 0,1 < j < n, according to 

the probability distribution 

(n - 1)! dpi... dp„_i. 

(c) Let a random vector x G be chosen according to the normalized uniform 

distribution the unit sphere {x G : ||a:|| = 1}. Then 

(M1) M2) • ■ • ) h"n-li h-m Mn-l-l) = (^^i -j- X 2 ,X^ + X^, . . . , X 2 n -3 + ^2n-2J ^2n-lJ ^2n) 
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are jointly distributed on the simplex h-j = 1 ) ^ Oj 1 ^ J ^ ''^+ !> aceording 

to the probability distribution 

(d) Let a random vector x G be chosen according to the normalized uniform 

distribution the unit sphere {x G : ||a;|| = 1}. Then 

{pi, ■ ■ ■ , Pm pn+l) {^1 ~t' 2 ^ 2 ) X^ + Xj^, ■ • ■ , ^2n—l ^2m ^2n+l) 

are jointly distributed on the simplex Ej=i Tj = 1 ; Mi ^ 0 ; 1 ^ i < 21 + 1 , according 
to the probability distribution 

r(«+2) 1 j j 

r(i) ^ji;^ dp,...dpm 

(e) Let a random vector x G be chosen according to the normalized uniform 

distribution the unit sphere {a; G : ||a;|| = 1}. Then 

( \/2i2i2i2 2 i2 i2 |2^ 

(M1) ■ ■ • ) Tnj (2^1 + 2;2 + 2:3 + 2:4,. . . , 2;4„_3 + 2:4„_2 + 2;4„_4 + 2:4„) 

are jointly distributed on the simplex Ei=i Mj = 1 ) Tj ^ 0,1 < j < n, according to 
the probability distribution 

n 

( 2 n - 1)! Pjdpi .. .dpn-i- 

i=i 

Proof. Part (b) has been proved in Cor A.3]. 

Part (a) follows from [22l Cor A.l] and a change of variables. The normalization 
constant comes from [22l Cor A.4]. 

Parts (c), (d), (e) can be proved using the arguments in the proof of [HI Cor A.3] 
and then using [HI Cor A.4] for the normalization constants. □ 


Appendix D. Jacobian Determinants 


Let us define the Jacobian determinants of real and complex maps. 

Recall that for a map / : M" —>■ K", (zi,..., z„) i-A (/i,..., /„), the Jacobian 
determinant is defined to be 


\d{fu- 

■Jn)' 

det 

'dfj' 

[a(zi,. 



_dzk_ 


To avoid confusion, for complex maps we will use a different notation. If each 
component of the map / : C" —?> C", (zi,..., z„) i-A (/i,..., /„) is analytic, then 
one can define the complex Jacobian determinant 


Jc 


\d{h,. 

■ ,/«)' 

det 

'df,' 

[a(zi,. 

■ 1 ^n)_ 


_dzk_ 


For our purposes however, the following will be more natural. For a (not nec¬ 
essarily analytic) map / : C" —>■ C", (zi,..., z„) i-A- (/i,..., /„), the real Jacobian 
determinant Jr is defined to be 



■Jn)' 

T 

’9(Re/i,Im/i,. 

. . , R,6 fri', IlR fn^ 

_d{zi ,.. 

■ 5 ^n). 


9(Rezi,Imzi,. 

.. ,Rez„,Im z„)_ 


We show below that for a given map /:€"—>■ C", (zi,..., z„) i-)- (/i,..., /„), 
its Jacobian Jr/ is equal to the Jacobian Jc/, where / : C^” —>■ C^" is the mapping 
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that takes (zi, zi,..., Zn, Zn) into (/i, /i,..., /„, /„), where we treat zj and zj as 
independent variables. 

Lemma D.l. For a map f : C” -)■ C”, (zi,...,z„) i-A 


d{zi,...,Zn) 


= Jc 


djfl, fl, ■ ■ ■ , fn, fn) 
d(^zi^ zi,..., Zfij Zji^ 


Remark. This implies that for an analytic map / : 


Proof. Indeed, 


C’’ 



■JnY 


7 

'difi,...,f„)) 

_a(zi,. 

■An). 



i9(zi,... ,z„)J 


dfj dfj 
dzk dzjk 

dzk dzk J l<j,k<r. 


= diag 


1/2 -i/2 

i/2 i/2 


d Re fj d Im fj 
d Re Zk 
dKcfj 


d Re Zk 
_ dlmj^ 

d Im Zk d Im Zk 


( diag 


l<j,k<n 


1 

—i 


where diag [A] denotes a block diagonal matrix with n copies of A on the diagonal. 

□ 


In the following lemma we let 


$(0) = Z” + KiZ^~^ + . . . + Kn-lZ + = ]^(z - Zj), 

4 = 1 

where we order z^-’s, e.g., by the absolute values and then by the argument. The 
mapping from the zeros to coefficients is then injective. It can be viewed as M" —>■ 
K” or as C" ^ C". 


Lemma D.2. 


(b) As a C" 


Proof. Let 


and 


(a) As an M" —>■ M" map, 

'd{Ki, ...,Kn) 


J 


9(zi,...,z„) 
C" map, 

'd{Ki,...,Kn) 


= A(zi,... ,z„); 


il(zi , . . . , Zyj) 


= |A(zi,...,z„)p. 


Sk= Y. 

ji<...<jk 


Pk=Y^3- 

i=i 

Then by the Newton identity, ksk = (—+ J2jZi{~^y~^Sk-jPj, and so 
inductively we get that ksk = (—l)^“^pfc+linear combination of si,..., Sfc-i- Thus 
noting that Kk = (—l)^Sfe and using Gaussian elimination, we obtain part (a) of 
the lemma. 

Part (b) then follows by using the remark after Lemma D.l and analyticity of 
the change of variables. □ 
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